REVISITING NUCLEOSYNTHESIS CONSTRAINTS 
ON PRIMORDIAL MAGNETIC FIELDS 



Dario Grassed and Hector R. Rubinstein 

Department of Theoretical Physics, Uppsala University 
Box 803 S-751 08 Uppsala, Sweden, and 
Department of Physics, University of Stockholm 
Vanadisvdgen 9, S-113 46 Stockholm, Sweden. 

Abstract 

In view of several conflicting results, we reanalyze the effects of mag- 
netic fields on the primordial nucleosynthesis. In the case the magnetic 
field is homogeneous over a horizon volume, we show that the main effects 
of the magnetic field are given by the contribution of its energy density 
to the Universe expansion rate and the effect of the field on the electrons 
quantum statistics. Although, in order to get an upper limit on the field 
strength, the weight of the former effect is numerically larger, the latter 
cannot be neglected. Including both effects in the PN code we get the 
upper limit B < 1 X I0 11 Gauss at the temperature T = 10 9 °K. We 
generalize the considerations to cases when instead the magnetic is inho- 
mogeneous on the horizon length. We show that in these cases only the 
effect of the magnetic field on the electrons statistics is relevant. If the 
coherence length of the magnetic field at the end of the PN is in the range 
10 <C Lq <C 10 11 cm our upper limit is B < 1 x 10 12 Gauss. 
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The study of the effects of magnetic fields on the primordial nucleosynthesis 
(PN) began with the pioneering works of Greenstein and Matese and O'Connell 
[0. In a previous paper Matese and O'Connell || first observed that magnetic 
fields larger than the critical value B c = eB /m e = 4.4 x 10 13 Gauss can affect the 
/3-decay rate, mainly through the modification of the phase space of the electrons 
and positrons. Greenstein advanced that, if magnetic fields were present in the 
early Universe, then this can affect not only the weak processes regulating the 
neutron to proton ratio in the early Universe, namely 

n + e + # p + 17 

n + v # p + e~ 

n # p + e~ + V , 

but also the expansion rate of the Universe. These effects will have consequences 
on the PN. However, they are competing effects. In fact, whereas the magnetic 
field energy density pb = B 2 /8tt accelerates the Universe expansion and therefore 
increases the predicted abundance of the relic 4 He, the effect of the field on the 
weak processes is to increase the n —>■ p conversion rate reducing the predicted 
4 He abundance. Greenstein argued that the former effect dominates over the 
latter. In a following paper, Matese and O'Connell || improved Greenstein's 
considerations by computing explicitly n — *> p conversion rate at the PN time in 
presence of a strong magnetic field and they confirmed Greenstein's qualitative 
conclusions. The comparison of the PN predictions for the relic light isotope 
abundances with the observations provides a criteria to bound the magnetic field 
strength at the nucleosynthesis time M. However, Matese and O'Connell did not 
provide any upper limit. 

Recently, the interest about magnetic fields in the early Universe received new 
impulse from the idea that these fields might be produced during the electroweak 
phase transition H or inflation |J. The renewed interest in this subject led several 
authors to investigate Matese and O'Connell's treatment again. 

Cheng, Schramm and Truran |7j first used the standard nucleosynthesis code 
[B] to get an upper limit on the magnetic field strength. The effects of the field 
they considered are the same studied by Matese and O'Connell. However Cheng 
et al. disagree with Matese and O'Connell in their conclusions. In fact, they claim 
that the effect of the magnetic field on the weak reactions is the most important. 
Cheng at al. obtained the following upper limit: B < 10 11 Gauss at the end of the 
nucleosynthesis (T = 10 9 °K). In both the Matese and O'Connell and Cheng at 
al. papers the magnetic field was assumed to be uniform over an horizon volume. 
We will return below to the crucial implications of this assumption. 

Recently, we reconsidered this subject ||. Besides the effects of the magnetic 
field on the reaction rates, we considered the effect of the field on the statistical 
distribution of the electrons and positrons. In fact, due to the effect of the mag- 
netic field on the electron wave function, the phase space of the lowest Landau 
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level is enhanced. If the temperature is small with respect to the energy gap 
between the lowest and first exited level, that is if eB <C T 2 [[ a relevant frac- 
tion of electron-positron pairs will condense in the lowest Landau level affecting 
the statistical distribution of the electron gas. As a consequence, the number 
density and energy density of electrons and positrons are increased with respect 
to the case where the magnetic field is not present. The expressions for the 
electron+positron energy density and pressure are 
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Here 7 = B/B c , e and are the electron, or positron, energy and chemical 
potential in mass units, and n labels the Landau level. 

It is evident that the modified contribution of electrons and positrons to the 
total energy density will affect the expansion rate of the Universe. The entropy 
content of the Universe will be also changed by the effect on the field on the 
electrons thermodynamics. However, for values of the ratio eB /T 2 ~ 1, the main 
effect of the modified electron statistics on PN is on the variation that it induces 
on the time derivative of the photon temperature. This derivative is given by 

dt dpem/dT 

where p em = p e + p 7 , p em = p e + p 7 , H is the Universe expansion rate and T is 
the photon temperature. For small values of the ratio eB/T 2 , the most relevant 
effect of the magnetic field enters in the derivative dp em / dT^ that is smaller than 
the value it would have if the field were not present. More physically, this effect 
can be interpreted as a delay in the electron-positron annihilation time induced 
by the magnetic field. This will give rise to a slower entropy transfer from the 
electron-positron pairs to the photons, then to a slower reheating of the heat bath. 
In fact, due to the enlarged phase-space of the lowest Landau level of electrons 
and positrons, the equilibrium of the process e + e~ <->• 7 is shifted towards its 
left side. As we already showed in ref.|J, and we are going to discuss in more 
details here, this effect has a clear signature on the Deuterium and 3 He predicted 
abundances besides that on the 4 He abundance. 



1 Observe that, neglecting dissipation, eB/T 2 remains constant during the Universe 
expansion. 
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Including the effect of the magnetic field on the electron statistics, together 
with the effect of the field energy density and the effect on the weak processes in 
the standard nucleosynthesis code we got the upper limit B < 3 x 10 10 Gauss. 
We also considered other kind of effects like that of the magnetic field on the 
nucleon |10[ and electron masses, but these effects proved to be subdominant. 

Recently our conclusions have been put in question in a preprint by Kernan, 
Starkman and Vachaspati [|llj. In this preprint, the effect of the magnetic field 
on the PN has been parameterized in terms of an effective number of neutrinos 
in order to avoid the numerical PN computations. Kernan et al. considered, at 
least in principle, all the effects of the magnetic field that we discussed above, 
including that on the electron statistics. However, although for opposite reasons, 
their main conclusion disagree both with Cheng et al. and our results. In fact, in 
Kernan et al. opinion the most relevant consequence of the magnetic field on the 
PN is the one mediated by the effect of the field energy density on the expansion 
rate of the Universe. The upper limit on the strength of the magnetic field they 
got is eBjTl < 1 that can be read B < 1 x 10 11 at T = 10 9 °K. This result does 
not differ from Cheng at al. numerical result although it does in its physical 
interpretation. 

The aim of this letter is to try to clarify this confused situation. As we are 
going to show, some mistakes are probably present in all refs. [[5J |T| || . To start 
with, we have recently realized that our upper limit, in ref.[[J is incorrect. This 
mistake was due to the use of a wrong version of the standard nucleosynthesis 
code. We have now rerun the right version of the PN code [F2 with an improved 
subroutine that evaluates the effect of the magnetic field on the electron and 
positron energy density, pressure and energy density time derivative. We consid- 
ered also the effects of the magnetic field on the weak reaction and the effect of 
the field energy density on the Universe expansion rate. Including all possible 
effects in the PN code we get the predictions we report in table I. 

These results have been obtained using N v = 3, the lowest value of the neutron 
half life compatible with experimental data, (Ti/2(n) = 885 s.) and the smallest 
value of the baryons to photons ratio that makes the predicted (D+ 3 He) /H ratio, 
in absence of the magnetic field, compatible with the astrophysical observations 
(rj = 2.8 x 10~ 10 ). This choice assured the minimal predicted abundance of the 
4 He for each value of the magnetic field strength that we considered. Requiring 
that such abundance do not exceed the value of 0.245 we infer from our data 
the upper limit 

B <lx 10 11 Gauss (4) 

when the temperature was T = 10 9 °K (end of PN) . 

In order to discriminate the relative weight of the effects we considered we 
have run our code leaving on only one effect at a time. We start considering only 
the effect of the magnetic field energy density. As expected, we see in Tab. II that 

2 However we used here an old version of the code, sec M. 
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Table I. Here we report the PN predictions for the relic light element abundance 
for several values of the parameter 7 = B / Be given at the temperature T = 
10 9 °K. 
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the predicted 4 He abundance grows with the field strength. In fact, larger values 
of B induce a larger Universe expansion rate then a higher freeze-out temperature 
for the n/p ratio. 

In Tab. Ill we report instead our predictions in the case only the effect of 
the magnetic field on the electron quantum statistic is considered. Note that, in 
order to display the relative weight of the effects that we considered, deriving the 
results of Tab. II and III we kept the values of 71/2(71) and rj fixed and equal to 
the value we used for Tab. I. 

It is interesting to observe the behavior of the sum of the predicted abundances 
for the Deuterium and 3 He while varying the magnetic field strength. This is 
qualitatively different from what reported in Tab. II. In fact, if only the effect on 
the electron statistics is considered, the (D + 3 He) /H abundances ratio decreases 
when the field strength increases (see also Fig.2). At the same time, the predicted 
4 He increases increasing B (see Fig.l). 

From the results we reported in Tab. II it follows that the upper limit we 
would get if only the effect of the magnetic field energy density were present is 
B < 2 x 10 11 Gauss whereas considering only the effect on the electrons quantum 
statistics (see Tab. Ill) we would get B < 5 x 10 11 0. Thus, contrary to our 
previous claim ||, the former effect numerically dominates the latter. However, 
from Fig. 1 it is clear at a glance that the effect of the magnetic field on the 
electron statistics cannot be neglected and indeed it dominates for small values 

3 We have to point out that this limit has been obtained using only our predictions for 
the 4 He abundance keeping the value of rj fixed. However, since 77 is an unknown parameter, 
it should be rescaled for every chosen value of B in such a way to get the largest D + 3 He 
abundance compatible with the observational data (see below). 
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Table II. Here we report the PN predictions obtained considering only the effect 
of the magnetic field energy density on the Universe expansion. 



7 (T 


= 


10 9 


°K) 


4 He 


(D+ 3 He)/H 


7 Li/H 











0.237 


1.04 x 10" 4 


1 


15 x 10" 


-10 


1 


X 


10" 


-3 


0.237 


1.05 x 10" 4 


1 


15 x 10" 


-10 


2 


X 


10" 


-3 


0.238 


1.06 x 10" 4 


1 


15 x 10" 


-10 


3 


X 


io- 


-3 


0.240 


1.09 x 10" 4 


1 


15 x 10" 


-10 


4 


X 


io- 


-3 


0.241 


1.12 x 10" 4 


1 


16 x 10" 


-10 


5 


X 


10" 


-3 


0.244 


1.17 x 10" 4 


1 


17 x 10" 


-10 


6 


X 


10" 


-3 


0.246 


1.22 x 10~ 4 


1 


19 x 10" 


-10 


7 


X 


10" 


-3 


0.249 


1.28 x 10" 4 


1 


23 x 10" 


-10 


8 


X 


10" 


-3 


0.252 


1.36 x 10" 4 


1 


28 x 10" 


-10 


9 


X 


10" 


-3 


0.256 


1.44 x 10" 4 


1 


34 x 10" 


-10 


1 


X 


io- 


-2 


0.259 


1.53 x 10" 4 


1 


42 x 10" 


-10 



of the field strength. In this sense we still disagree with the conclusion of the 



authors of ref.flll]]. It may be that in ref. [0] the effect of the magnetic field on 
the electrons and positrons statistical distribution was not treated properly close 
to the saturation of the upper limit on B, that is for eB/T 2 ~ 1. 

Concerning the effect of the magnetic field on the weak reactions we have 
verified that such effect is negligible with respect to the others. Mass changes are 
also negligible in the range of field strengths that we considered. About these 
conclusions we agree with the authors of the refs. M and [FT] . Since the numerical 
predictions of ref. do not differ significantly from those of ref. |TT| and the results 
we reported in Tab. II we think that the different physical interpretation of the 
results given in ref.JFj is probably due to an oversight writing that paper. 

To summarize, we have found that the main effects of a magnetic field on the 
PN plays through the action of the field on the electrons and positrons quantum 
statistics and the direct effect of the field energy density on the expansion rate 
of the Universe. Although the latter effect is numerically larger than the former, 
these effects are indeed comparable, at least if the magnetic field is uniform over 
the horizon scale as we assumed so far. 

Let us now discuss how our previous consideration will be affected if we relax 
the uniformity assumption. Before discussing more model dependent motivations 
to consider inhomogeneous magnetic fields, it is mandatory to compare the upper 
limit we got above with the upper limit one would get considering the effect of 
an homogeneous cosmic magnetic field on the early Universe isotropy. In fact, 
since magnetic fields break rotational invariance they can induce an anisotropy 
in the Universe expansion. This anisotropy can have observable consequences on 
the isotropy of the cosmic background radiation and on the PN [14}|[115"|. Using 
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Table III. Here we report the PN predictions obtained considering only the effect 
of the magnetic field energy density on the electrons and positrons quantum 
statistics. 
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these criteria Zeldovich and Novikov [|14[] got the limit B < 4 x 10 Gauss for 
a magnetic field homogeneous at least over a horizon volume considered at the 
temperature T = 10 9 K. This limit is roughly one order of magnitude more 
stringent than the one we got above as those obtained in refs.|7[] and (ll[as well. 

Nevertheless, if the magnetic field is inhomogeneous with a coherence length 
much smaller than the Hubble radius we do not have to care about the Zeldovich 
and Novikov upper limit. In our opinion this is the most plausible physical sce- 
nario. It is in agreement with the predictions of most of the models for magnetic 
field generation in the early Universe [plj. Among these models we are mainly 



interested in those that predict that magnetic fields are generated during the 
electroweak phase transition. The predicted coherence length of the field at the 
electroweak phase transition time is L Q rs m^ 1 ~ lO -16 ^^ 1 (where is the 
Hubble radius at that time) in the case the transition is second order || or, per- 
haps a more probable scenario, L ~ Rbubbie ~ ^ ^^H^ for a first order phase 
transition |r7J[TJ|]. Assuming the magnetic field strength decreases only due the 



Universe expansion the ratio Lq/H -1 remains constant in time. In both cases it 
is evident that L H* 1 . 

Having magnetic fields that fluctuate in space on a typical scale Lq H^ 1 
forces us to reconsider the analysis we made in the first part of this letter paying 
some attention to the different scales on which the effects of the field on the PN 
act and the nature of the averaging processes. 

To start with, let consider the effect of the field on the weak processes. The 
mean free path between two weak reactions is Tw ~ G 2 F T 5 that is of the order 
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Figure 1: The 4 He predicted abundance is represented in function of the param- 
eter 7, considered at T = 10 9 °K, in three different cases: only the effect of the 
magnetic field energy density is considered (dashed line); only the effect of the 
field on the electron statistics is considered (dotted-dashed line); both effects are 
considered (continuous line). The dotted line represents the observational upper 
limit. 



of the horizon radius at the onset of the PN {H~ X {T ~ 10 10 °K) ~ 10 10 cm). It 
is clear that if L <C T^ 1 the charged particles involved in the process will feel 
a magnetic field changing a number of times between the two reactions so that 
the effect of the field will be almost averaged out [fT9[] . At any rate, the effect of 
the field on the weak reaction rates was already subdominant in the case B were 
assumed to be homogeneous. 

More relevant is, instead, the effect of the inhomogeneity of the magnetic field 
on the action that the magnetic field energy density p% plays on the Universe 
expansion. As discussed in detail in refs. [^9[pD| , since in this case the field is 
fluctuating in space, only a mean magnetic field within a horizon volume can be 
defined as a meaningful quantity. This quantity is defined by 

Pb = T^T d 3 rBl ls (r) (5) 



2V H Jl 

where 



B m . = B (£.) (^) (6) 



is the root-mean-squared field strength. Here p is a parameter that depends 
on the statistical configuration of the magnetic field. Typical values of p are 
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Figure 2: The (D+ 3 He)/H predicted abundance is represented for the same 
cases illustrated in Fig.l. 

1/2, 1, 3/2. Inserting Eq. @ in Eq. (H) it is evident that ps will suffer a huge 
suppression if L <C so that it will have no chances to affect any more the PN. 
The impact of the magnetic field on the electron quantum statistics is much less 
affected than the others by the inhomogeneity of the field. In fact, in this case, the 
characteristic length scale we have to compare to Lq is the Compton scattering 
length of the electrons Ac*, since this scattering is responsible for the electron 
thermalization, and the radius of periodic motion of the electrons in the plane 
normal to the magnetic field vector. At the PN temperature (T ~ 10 9 ^ 10 °K) 
the Compton cross section does not differ significantly from the Thomson cross 
section. Then we have 

A c = (a T n e v)- 1 ~ 1(T (2 ^ 5) cm. (7) 

This has to be compared with the Hubble radius at the same temperature that is 
~ lO 10 ^ 11 cm. The radius of the classical orbit of the electron in an over critical 
magnetic fields is given by 

^n=(^) 1/2 m e - 1 (8) 

where n labels the Landau level. Since the magnetic field plays its main effect 
of on the lowest Landau level it is clear that the radius we have to care about is 
much smaller than Ac and far below the expected value of Lq at the PN time. 
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Since, Ac <C T^ 1 <C H^ 1 we see that the effect of the magnetic field on the 
electron quantum statistics provides the best probe of the field at small length 
scales. If Ac <C Lq <C H^ 1 this is the only effect we remain with. 

The PN predictions of Tab. Ill cannot be used in their present form to deter- 
mine, in this case, the upper limit on the magnetic field strength. In fact, as we 
pointed-out above, obtaining Tab. Ill we kept the value of rj fixed. However, 77 
is an unknown parameter and its value has to be bound from below in order the 
predicted abundances of D and 3 He do not exceed the relic abundances extrap- 
olated from present observations. Since in our case the predicted of abundances 
of these isotopes depend on the field strength, we have to rescale 77 for every 
value of B in order to get their largest predicted abundances compatible with 
the observational upper limit. According to ref. we assume this limit to be 
(D+ 3 He)/H < 1.1 x 10 -4 . Finally, we have to compare the predicted abundance 
of the 4 He, obtained using such value of 77, with the maximal abundance compat- 
ible with observations 4 He < 0.245 |i~3| . Using these procedure we get the upper 



limit 

B < 1 x 10 12 Gauss. (9) 

Again, this limit refers to the temperature T = 10 9 °K. The corresponding 
value of rj is 2.3. As expected, this value is smaller than the standard lower 
limit rj > 2.8, since the effect of the magnetic field on the electron statistics 
suppresses the predicted abundance of D and 3 He. It is worthwhile to note 
here that, although this effect increases also the predicted abundance of 4 He, 
large magnetic fields might help to alleviate the "nucleosynthesis crisis" that is 
presently under debate [^1[ . 



So far we neglected any dissipative effect in the plasma. However, even in a 
relativistic plasma conductivity a is a finite quantity. This means that any field 
configuration on a length scale L <C Ld iss , where 
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will be dissipated. It is interesting to observe that if eB /T 2 ^> 1 the conductivity 
properties of the plasma will be seriously affected by the magnetic field. In fact, 
in this case the e + e~ pair number density is increased by the effect of the field 
on the electrons and positrons phase-space. As a consequence this will make the 



plasma conductivity to decrease then L<a ss to grow []25 . This can have important 
consequences on the field evolution after it is generated in the early Universe. 
However, such effect is negligible at the PN time. Following ref. |20| Ldi SS at 
T = 10 9 °K MeV can be estimated to be 
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0.1^ 1/4 — « 1 10 cm. (11) 
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Taking this considerations into account we conclude that our upper limit (|9|) 
applies only if 10 <C L <C 10 11 cm. 

It is worthwhile to see how our upper limit constraints at least one of the 
recent models of relic magnetic field generation. Assuming the field is generated 
at the end of a first order electroweak phase-transition [r?| ||18|| , within large the- 
oretical uncertainties the predicted magnetic field strength at the temperature 
T ~ 10 15 °K is « 10 24 Gauss with a coherence length L m 10~( 2 ^ cm. If the 
magnetic field flux above this length scale remains frozen in the plasma, so that 
the field strength decrease only due to the Universe expansion, at the PN time 
we expect: BiT = 10 9 °K) ~ 10 12 Gauss and L ~ 10 4 ^ 5 cm. Thus, this model 
is not excluded by PN considerations. 

In conclusion, in this letter we unfortunaly revise upwards previous limits 
on homogeneous magnetic field strengths at PN time by a factor 10. Earlier 
calculations contained mistakes on a variety of points as discussed. As a conse- 



quence, previous bounds based on Universe isotropy remain more stringent 1 14 



However, in the case the magnetic field is inhomogeneous on the horizon length 
scale, isotropy considerations do not apply. In this case, we showed that the main 
observable effect of the magnetic field on the PN is on the electron and positron 
quantum statistics. 
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